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Abstract 

In our previous papers we proposed a continuum model for the 
dynamics of the systems of self-propelling particles with conservative 
kinematic constraints on the velocities. We have determined a class of 
stationary solutions of this hydrodynamic model and have shown that 
two types of stationary flow, linear and radially symmetric (vortical) 
flow, are possible. In this paper we consider the stability properties of 
these stationary flows. We show, using a linear stability analysis, that 
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the linear solutions are neutrally stable with respect to the imposed 
velocity and density perturbations. A similar analysis of the stability 
of the vortical solution is found to be not conclusive. 

1 Introduction 

The dynamics of the systems of self-propelling particles (SPP) is of a great 
interest for physicists as well as for biologists because of the complex and 
fascinating phenomenon of the emergence of the ordered motion. In nature 
these systems are represented by flocks of birds, schools of fishes, groups of 
bacteria etc. pQ-[3]. From the physical point of view many aspects of the 
observed non-equilibrium phase transition from disordered to ordered motion 
are to a large extent still an open problem. 

The first numerical model simulating the behavior of the SPP was pro- 
posed by T. Vicsek at al. [4J. Their model is based on a kinematic rule 
imposed on the orientations of the velocities of the self-propelling particles. 
At a low noise amplitude and a high density it was shown that the system 
undergoes the transition from the disordered state to coherent motion. The 
nature of the transition is not established yet. T. Vicsek's investigation shows 
that the occurring transition is of the second order. Several extensions of T. 
Vicsek's model have been proposed, which consider particles with varying ve- 
locities, different types of noise, and including external force fields and/or in- 
terparticle attractive and repulsive forces [5]- [7]. The simulations performed 
in [7J show the discontinuous nature of the transition in the T. Vicsek's model. 

Properties of T. Vicsek's model were also investigated analytically. In [5] 
the spontaneous emergence of ordered motion has been studied in terms of 
so-called control laws, using graph theory. Generalizations of the control laws 
were considered in [91 [TUJ. In [10J it was shown that the organized motion 
of SPP with the control laws depending on the relative orientations of the 
velocities and relative spacing, can be of two types only: parallel and circular 
motion. The stability properties of these discrete updating rules (including 
the T. Vicsek's model) and the dynamics they describe were considered using 
Lyapunov theory in [HI [9], [11] . 

In our first paper [12J we constructed a hydrodynamic model for the 
system of self-propelling particles with conservative kinematic constraints, 
which can be considered as a continuum analogue of the discrete dynamic 
automaton proposed by T. Vicsek et al. 
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Based on the conservation of the kinetic energy and the number of parti- 
cles our model is represented by the following equations: 



= u,(r,t)xv(r,t), (1) 
^ + V-(n(r,f)v(r,t)) = 0, (2) 

where v (r, t) and n (r, t) are the velocity and the density fields respectively 
and uj (r, t) is an angular velocity field which takes into account the non- 
potential character of the interactions between the particles. We modeled 
this field as follows: 

u) (r, t) = J Kx (r - r') n (r', t) rot v (r', t) dr' + 

J K 2 {r- r') Vn (V, t)xv (r', t) dr', (3) 

where Ki 2 (r — r') are the averaging kernels. In particular we considered a 
simple case of averaging kernels: 

Ki (r — r') = Si 5 (r — r') , where z = 1 or 2. (4) 

We call this the local hydrodynamic model (LHM). In this case Eq. ([3]) 
reduces to 

uj (r, t) — sin (r, t) rot v (r, t) + s 2 Vn (r, t) x v (r, t) , (5) 

where 

r Ki^dr. (6) 



In our second article [13] we have shown that the only regimes of the station- 
ary planar motion in our model are either of translational or axial symmetry. 
In this respect our continuum model gives results similar with those obtained 
in the discrete model of T.Vicsek [HE]. 

In this paper we investigate the stability of the obtained regimes of motion 
with respect to small perturbations. In the next section we consider the 
stability of the planar stationary linear flow with respect to the velocity 
perturbation directed along the stationary flow and perpendicular to the 
flow. We show that in both cases the evolution of the perturbations has 
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an oscillatory behavior, which means that they neither grow nor decay with 
time. This can be interpreted as neutral stability [H] of the corresponding 
stationary flow. Also the external pressure term —Vp/n can be included into 
Eq. ([I]) in order to account for potential external forces. In such with 
S2 = there exists the special case of the incompressible flows, n = const, 
when the equations of motion coincide with that for potential flow of ideal 
fluids. As is known [H] such motion in 2D geometry is stable in the Lyapunov 
sense. 

In the third section we consider the stability of the planar stationary 
radially symmetric (vortical) motion of SPP with constant velocity and the 
density. We find that in this case the linear analysis does not lead to a 
conclusive answer about the stability of the solution. 

2 Stability of planar stationary linear flow in 
the local hydrodynamic model 

2.1 Stability with respect to a velocity perturbation 
along the flow 

In this section we consider the stability properties of planar stationary linear 
flow for the local hydrodynamic model with s 2 = 0, which we further call 
local hydrodynamic model 1 (LHM1). At the end of the section we will 
shortly discuss how these results extend to the local hydrodynamic models 
with si = and Si = s 2 - For LHM1 the stationary linear flow is given by 

v (r) = v e x and n (r) = n , (7) 

where vq and no are constants. 

We consider velocity and the density perturbations of the following form: 

vi(r,t) = v A||e ik - r e Q ii t e a . and n x (r, t) = n B\\ e ikr e""*, (8) 

The velocity perturbation chosen is directed along the stationary linear flow. 
Here Am, Bu are constants, k = k x e x + k y e y is the wave vector and an is an 
exponent, which determines the time evolution of the perturbation. 

Substituting the solution v (r, t) = v + Vi (r, t) , n (r, t) = n + n\ (r, t) 
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into Eqs. ([I])-© we obtain the linearized system of equations: 

<9vi 



at 



+ (v -V)vi = sirio (rot vi) x v , (9) 



+ V ■ (novi) + V • (mv ) = 0. (10) 
For the perturbation (|SJ) this system reduces to 



<9y 



0, (12) 



-sr + ^— + ^o^- = 0. (13) 
at ox ox 

Using Eq. (JHJ) one may obtain the relation between an and the wave number. 
From Eq. (Ilip it follows that 

«H = -ik x v (14) 
whereas from the linearized continuity equation (fl3|) we have 



ot\\ = -ik x v 



Both the equalities are satisfied only in the case when A\\ = 0. 

Thus, in the linear stability analysis with respect to small deviations of 
the velocity and density fields, we obtain the following perturbed solution 

v = v e x , n = n [l + B\\ e lk * y e ik ^ x ~ Vot) ] . (15) 

Taking the real part of the density perturbation we have 

v = vq e x , n = n [l + B\\ cos (k • r — k x Vot)\ . (16) 

The corresponding density field is shown on Figl. 

This flow Eq. ffl6|) should satisfy the linearized system of the constraints 
(conservation of the kinetic energy and the number of particles) which are 
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Figure 1: The total density field n(r,t) /n and the stationary solution 
n/riQ = 1 as a function of x* = k x x and t* = k x Vot for k y = 0. 

imposed on any solution of our model. This implies that the following con- 
ditions must be fulfilled: 

J m (r,t)dr = 0, (17) 
J [2n (v -v 1 (r,t))+n 1 (r,t)v 2 ]dr = 0. (18) 

Since Vi (r, t) — both conditions reduce to 

/ rei (M)*^ B|l /e^/e'M^» (i , = . (19) 

If one integrates Eq. [19] over the period of the integrand, one may can that 
this condition is fulfilled. 

The obtained perturbed flow is an oscillatory field (perturbation oscillates 
with a frequency a\\ as t — > oo) which means that the corresponding station- 
ary solution is neither stable nor unstable within the first order perturbation 
theory. In other words we may conclude that in our local hydrodynamic 
model the stationary linear flow is neutrally stable with respect to a small 
density field perturbations. 

The stability analysis of the other possible hydrodynamic models with 
si = or s\ = S2 gives qualitatively similar result. 



6 



2.2 Stability with respect to a velocity perturbation 
perpendicular to the flow 

In this section we investigate the stability properties of the stationary linear 
flow in the LHM1, Eq. (j7|), with respect to a velocity perturbation normal 
to the stationary flow. We consider only a velocity perturbation, which we 
take in the form of a plane wave: 

V! = v A± e ik r e a±t e y} n x = 0, (20) 

where Aj_ is a constant, k is a wave vector and the exponent a± describes 
the time evolution of the perturbation. 

Substituting the perturbation in the linearised equations (Ml- CTUl) it fol- 
lows that 

dvi , , dv\ . . 

-gl+Ml-in.)-^ = 0, (21) 

= and L = 0, (22) 

dy 

which imply that 

a ± = ik x v (si7i - 1) . (23) 

Thus the time evolution of the perturbed velocity field is determined by the 
purely imaginary exponent in Eq. (|23|) : 

v = v + Vl (x, t) = v [e x + A ± e tk ^ x+V ^e y ] , n = n (24) 

where the "phase speed" is given by 

V = v (sin - 1) . 

Taking the real part of the velocity perturbation we obtain as final result: 

v = v [e x + A ± cos [k x (x + Vt)} e y ] , n = n . (25) 

The corresponding velocity profile is shown in Fig.2. 

Since the velocity perturbation was taken to be normal to the unperturbed 
field and n\ = 0, both of the constraints of the constancy of the kinetic energy 
and the number of particles, Eqs. (TTTj) - (TTH]) . are satisfied. 
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Figure 2: The total velocity field v (x, t) /t> and the stationary velocity field 
Vo/vq = 1 as a function of x* = k x x and t* = k x Vt. 

As one may see the time dependent part of the velocity perturbation is 
a finite oscillatory function which means that the corresponding stationary 
solution is neutrally stable. 

As in the previous section the stability analysis of the other possible 
hydrodynamic models with s\ = or s\ = S2 gives qualitatively similar 
result. 

3 Stability of stationary vortical flow with 
constant velocity and density in the local 
hydrodynamic model 

As we have shown in our previous article [13j there are two classes of the 
stationary flows in the LHM, linear and radially symmetric or vortical flow. 

The stationary vortical solution of the LHM1 (s2 = 0) is given by vo (r) = 
v v (r) e v , n Q (r) = n (r), [12], where 



Here ro is a cut-off radius of the vortex core and the constant C s t is detemined 




(26) 
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by the circulation of the core 



vdi = a 



st ■ 



(27) 



r=ro 



We consider small perturbations Vi (r, if, t) of the velocity field and n\ (r, if, t) 
of the density field. The linearized system in the LHM1 is then given by 

dvi 



dt 
diii 



f (Vi • V) V + (V • V) V! 

+ V ■ (novi) + V ■ (nivo) 



si^o [(rot vi) x v + (rotvo) X vi] 
+sini (rotv ) x v , (28) 
0. (29) 



In this section we consider the stability of a particular class of station- 
ary vortical flow for which the density is constant and given by Uq = 1/si. 
Substitution in Eq. (|26|) results in a constant velocity field v = v^b v = 
(C s t/27iro) e v = Ce v . We write the small perturbation in the general form 

V! = a (r, <f, t) e r + b (r, if, t) e v and n\ = n Ci (r, if, t) . (30) 

For the projections of the velocity field v = vo (r) + vi (r, if, t) together with 
the continuity equation for the density field n = no + n\ (r, if, t) we have 

drb da 



da ^bv 
dt r 



v v da 
r df 
db v v db 
di^ 



dr dip 



i 



dci I 
dt r 



dra db 
dr dip 



+ 



r dip 
v v dci 
r dip 



0. 
0. 



- V~ Cl f ' (31) 



(32) 
(33) 



In order to simplify the problem we restrict our discussion to the case with the 
radial component of the velocity perturbation being constant, i.e. a (r, if, t) = 
const. 

Then one can transform equations (l3~T i) -( l33l) into 

db v v db 

dt r dip 
db Civ v 
dr 



0. 



dci | 1 
dt r 



a + 



db\ 



+ 



r 

v v dci 



dip J r df 



0. 



(34) 
(35) 
(36) 
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The velocity perturbation must be a periodic functions of the angle p and 
can therefore be written as: 



b(r,(p,t) =v ip B (r) e imp e^\ 



(37) 



where B (r) is a function of r, m is an integer and f3 is a constant factor, which 
describes the time evolution of the perturbation, Eq. fl30l . Substituting this 
into Eq. fl34|) one obtains 

(38) 
(39) 



a ■ v v 
p = —im — 

r 



and consequently 

b (r, (p, t) = v v B (r) exp 
From Eq. (|35|) it follows that 



im I (p t 

r 



Ci (r, v?,t) = -r 



OB (r) 
Or 



+ zto- 



( cp -t 



(40) 



Substituting this into Eq. ( |36l) we obtain that a (r, ip, t) = 0. 

The solutions ( 1391) and ( HU1) satisfy the linearized system of constraints, 
Eqs. ( 1T7|) and ( |T8|) . as one can see by angular integration. 
Thus, we see that the time evolution of the perturbation Eq. (!30|) is deter- 
mined by the purely imaginary exponent Eq. ( 1381) . 
Taking the real part in Eqs. (|39|) and (j4"0"]) we obtain 



b{r,(p,t) 
m (r,v?,t) 



(r) cos 



r 



n 



t sin 



m\p> 1 



dB(r) 

—r — cos 

or 



m I <p 



■t 



(41) 



(42) 



As a result the whole solution for the velocity and the density profiles has 
the following form: 



n(r,(p,t) 



1 + B (r) cos 



m l <p 



r 



n Q <l + 



mv,„ B (r) 



t sin 



m [p 1 

r 



OB (r) 
—r — cos 

or 



m [p —t 

< - r 



(43) 



(44) 
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Figure 3: The total velocity field v (r, <p, t) / v 9 and v v (r) /v^— 1 as a function 
of ip and t* = v v t / r for m = 1 and r = 5 m. 

The velocity field is shown in Fig. 3 for m = 1 and r = 5 m. 

Together with the oscillatory contributions we now also have the contri- 
bution proportional to t times an oscillating function. This does not neces- 
sarily mean that the stationary vortical flow is unstable. The obtained result 
indicates that the linear analysis does not give a definite answer regarding 
the stability of the stationary flow. The linear analysis does not give the 
definitive answer regarding the stability of the stationary flow and further 
investigation of higher order terms is required. This is beyond the scope of 
the present paper. 

4 Conclusions 

In this paper we considered the stability properties of the planar stationary 
flows of the local hydrodynamic model constructed in our first paper for 
a system of self-propelling particles [12J. These flows are the linear flow 
and the radially symmetric flow. Our analysis shows for linear flow, using 
linear perturbation theory, that the time evolution of the imposed velocity 
and density perturbations are oscillatory. It follows that the linear flows are 
neutrally stable. For radially symmetric (vortical) flow linear perturbation 
theory does not lead to a conclusive result. A definitive answer about the 
nature of the stability can only be given by considering also higher order 
terms in the perturbation expansion. Such an analysis is beyond the scope 
of the present paper. Note that such a situation is typical for Hamiltonian 
systems which are conservative by definition and therefore do not display an 
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asymptotic type of stability [13] . 
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